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Abstract 

The photon energy distribution in the inclusive 6 — s + 7 transitions is a combina- 
tion of two components: the first component, soft physics, is determined by the so 
called primordial distribution function, while the second component, perturbative 
physics, is governed by the hard gluon emission. A simple ansatz is suggested for 
the primordial distribution function which obeys the QCD constraints known so 
far. We then discuss in detail how the hard gluon emission affects the energy dis- 
tribution. An extension of the Sudakov approximation is worked out incorporating 
the Brodsky-Lepage-Mackenzie prescription and its generalizations. We explicitly 
calculate the marriage of nonperturbative with perturbative effects in the way re- 
quired by OPE, introducing separation scale fi. A few parameters, such as rub and 
/i^ affect the shape of the distribution and, thus, can be determined by matching 
to the experimental data. The data, still scarce, while not giving precise values for 
these parameters, yield consistency with theory: the current values of the above 
parameters lie within experimental uncertainty. On the theoretical side we outline 
a method allowing one to go beyond the practical version of OPE. 



1 Introduction 



The inclusive 6^5 + 7 transition is one of a few low-energy processes which is 
potentially sensitive to new physics. Although the recent observation of this transi- 
tion has so far detected no deviations from the Standard Model expectations, it 
is not ruled out that contributions to the decay parameters going beyond the Stan- 
dard Model will be discovered in the future 0. Therefore, untangling the effects 
of the electroweak interactions from those of the strong interactions in this problem 
becomes an important and urgent task. 

The purpose of this work is to develop a theory which is as close to fundamental 
QCD as possible. Some model dependence is, unfortunately, unavoidable in the 
photon energy distribution at the present level of the strong interaction theory. Our 
approach in this aspect is to provide guidelines that can readily accomodate future 
refinements, if necessary, rather than giving a final prescription. In other words, we 
try to set a general framework for the theoretical analysis of 6 — 57 and similar 
transitions. 

The basis of our approach is the heavy quark expansion, a version of the Op- 
erator Product Expansion (OPE)0]. This method has experienced a dramatic de- 
velopment over the last decade (see Ref. for a review and a list of references). 
The total inclusive widths of heavy flavor hadrons are most directly amenable to 
calculations within the 1 / mq expansion 0, |^ . In particular, it was shown that the 
nonperturbative corrections to the total widths start from terms ~ I/^^q [0- 
the radiative decay 6^5 + 7 they were first calculated in Ref. and turn out 
to be rather small. The absence of the l/rriQ corrections does not apply, however, 
to the shape of the photon spectrum which depends in a crucial way on details of 
the bound state, viz., on the nonrelativistic "Fermi motion" of the heavy quark in- 
side the decaying heavy hadron 0. This nonrelativistic 'jiggling' leads to essential 
nonperturbative effects in the kinematical region near the endpoint, of the relative 
width l/rriQ, in accord with the naive expectations following from the nonrelativistic 
picture. However, the way this "Fermi motion" manifests itself in QCD is essentially 
different ^ from the description appearing in the nonrelativistic picture. A gen- 
eral formalism appropriate for treatment of the endpoint domain of the spectrum in 
QCD was suggested in Refs. |jlO|, |ll|, |12|. This formalism was applied to the analysis 
of nonperturbative effects shaping the endpoint photon spectrum in the transition 
6 — s> 57 in Refs. [|^, 0, |I4|. The results can be summarised as follows: the spectrum 
is expressible in terms of a universal (one-dimensional) distribution function whose 
moments are related to local heavy quark operators of increasing dimensions. This 
distribution function replaces the Fermi motion wavefunction of the nonrelativistic 
quark model. In contrast with the latter, the properties of the distribution function 
crucially depend on whether the final state quark in the process at hand is light 
or heavy |jl2[. If the s quark mass is neglected the distribution function we deal 
with in 6 — > 57 is the light-cone one. The physical photon spectrum is obtained 
as a convolution of the primordial (soft) distribution function with the perturbative 
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spectral function corresponding to the hard gluon emission. 

Attempts to account for the motion of the heavy quark inside heavy hadrons 
at a purely phenomenological level were made previously (e.g. WB, ITBl). Models of 



this type [17] treat the heavy hadron as a bound state of the heavy quark Q plus 
a spectator, with a certain momentum distribution. Although at first sight this 
approach, heavily relying on the nonrelativistic picture, contradicts QCD, a closer 
inspection shows |jl9| (see also ||20|) that the models of Refs. [jl5|, [T^ are compatible 
with QCD through the leading order in l/rriQ for the decays into light quarks. 
The corresponding ansatz reduces to a specific choice of the primordial distribution 
function. In particular, the Gaussian distribution over the Fermi momentum in 
conjunction with the energy-momentum conservation relation of the quasifree type 
(routinely postulated in phenomenological models) implies the so-called "Roman" 



function for the primordial distribution |jT9| 



$(a;) = 1^(1 -x)e~H^^")-T^)' ; (1) 

the dimensioneless parameters a and jS are certain functions of Tn1^/p\ where mgp 
and pf are parameters of the AC^M^ model |jl6[. 



The models of Refs. [0, |l6l, being compatible with QCD in principle, are still 
unsatisfactory for a number of reasons. First, the Roman distribution function, 
taken at its face value, may have problems accomodating the values of the basic 
QCD parameters A and /i^, as they emerge from the recent estimates (see below). 
Second, these models per se are responsible only for a soft smearing of the photon 
spectrum. Effects due to the hard gluon emission in the decay process that smear 
the spectrum further have to be incorporated additionally. The proper combina- 
tion of the soft and hard effects in the spectrum crucially depends on the way one 
introduces the normalization point /i, an aspect which is usually totally ignored in 
the corresponding discussions. Finally, some technical details in the application of 
the models to analysing experimental data were not completely consistent with the 
heavy quark expansion (see Ref. for details). In this work we address these 
issues in turn. 

We propose a simple two parameter ansatz for the distribution function consistent 
with all general properties expected for the light-cone distribution function (which 
is relevant when the heavy quark decays into a massless one). Our ansatz seems 
to be more flexible allowing one to accommodate realistic values of the basic QCD 
parameters. We then dwell on the main aspects of inclusion of the hard gluon 
corrections. The endpoint behavior of the spectrum is affected by the perturbative 
radiative corrections which have a double logaritmic enhancement in the case of the 
light flnal quark (i.e. in the 6 — >■ S7 transition). Although for the actual b quark 
mass the logarithm of the mass is not a particularly large number, the presence of 
Sudakov's double logarithms enhances the role of the perturbative gluons. 

An interplay between the perturbative and nonperturbative smearing of the spec- 
trum is quite peculiar. A version of the Sudakov form factor appears that essentially 
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reduces the absolute height of the spectrum in the endpoint region. At the same 
time, the shape of the spectrum is formed mainly by the primordial distribution 
function [|T^]. Accounting for the bremsstrahlung corrections is crucial for centering 
the predictions around the realistic values of the hadronic parameters. 

Decreasing the normalization point /j, we change the definition of which gluons are 
to be treated as soft and which as hard. Accordingly, the primordial distribution and 
the hard perturbative one vary simultaneously - a part of the primordial smearing 
leaks into the perturbative smearing. In principle, the physical distribution function 
must be /i independent, of course. However, since approximations are made in 
different ways in the soft and hard domain, the n independence of the physical 
distribution function may hold only approximately. One of our tasks is to check this 
approximate /i independence. To achieve this goal we, obviously, must explicitly 
introduce the separation scale /i in calculating the perturbative distribution (i.e. 
the coefficient functions in the Wilson OPE). Thus, we go beyond the practical 



version of OPE 21 



After the theoretical framework is set we then compare the theoretical spectrum 
with the CLEO data The data, however, are extremely crude at present. Within 
the existing uncertainty the parameters cannot be precisely attained. Still, though, 
one can make some attempts with the data and show that the current level of 
analysis is not only theoretically consistent, but also consistent with information on 
the basic QCD parameters obtained from other sources. It then remains for the 
data to improve for the analysis to become more predictive. 

The organization of the paper is as follows. In Sect. 2 we formulate the problem 
and remind the basic elements of the theoretical description. In Sect. 3 an ansatz 
for the primordial light-cone distribution function is suggested. Sect. 4 is devoted to 
perturbative corrections; here we present both rather standard expressions and new 
considerations called upon to get a more accurate description in a manner consistent 
with OPE. In particular, in Sect. 4.6 a method is outlined allowing one to explicitly 
introduce the separation scale /x. Sect. 5 finalizes our expressions for the total 
spectrum. In Sect. 6 we briefiy discuss numerical results and give a comparison 
with the first experimenatal data available now. In Sect. 7 our conclusions are 
summarized. 



2 Formulating the Problem 

The 6 — > S7 transition at the fundamental level is generated by electroweak penguins 



2^ . Since in this work we are interested only in the strong interaction effects at 
distances m'^^ and larger we will just assume that the effective Lagrangian governing 
this transition has a local form 

A^s7 = ^F^^s{l + ^^)ia^,yh (2) 
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where his a constant which includes the electroweak physics as well as the hard gluon 
corrections p3[, both coming from distances smaller than m^J"^ [^. We hasten to 



add that the assumption of locality of the interaction (|^) is not absolutely correct. 
We will return to this point at the end of the Section. We also note that the s quark 
(current) mass will be consistently neglected. 
The total inclusive decay width is given by 



r(6 ^ S7) = To 

where 



1 + 71- + 72 (-)' + ■■■ + Oil/ml 



(3) 



7i,2,... are numerical coefficients, and all quantities, including h and m^, are assumed 
to be normalized at the scale mi,. The expansion in the square brackets is well-defined 
within Wilson's OPE. In particular, all coefficients in the perturbative expansion are 
infrared finite, and the whole series is well-behaved. We will not be interested in the 
corrections to the total width which, as can be explicitly seen, are small. 

The subject of our study is the photon energy spectrum which, unlike the total 
width, is strongly modified by corrections already at the level C(l/mQ). Indeed, in 
the parton model (valid in the academic limit mg oo) the adequate description 
of the inclusive transition is given by the two-body decay: a free b quark decays into 
a free s quark and a photon. Accordingly, the photon spectrum is monochromatic 

6(e-^) (4) 



TodE V 2 

where E is the photon energy. It is clear that in the actual decay of the B meson 
the monochromatic line aX E = mb/2 is smeared, and the photon spectrum stretches 
both, upwards, to the kinematic boundary at £^ = Mb/2, and downwards, to zero. 

Two distinct effect are responsible for this smearing. First, the decaying heavy 
quark is not at rest in the restframe of i?. It is submerged in a soft medium, light 
cloud, which exchanges energy and momentum of order of Aqcd with the heavy 
quark at hand. Thus, the soft gluon medium from the cloud creates a primordial 
nonperturbative distribution which is solely responsible, in particular, for filling in 
the window between mb/2 and Mb/ 2 |^. Second, in the process of decay the quarks 



can emit and absorb hard gluons, i.e. those whose momenta are much larger than 
Aqcd- For example, the s quark can shake off some of its energy and momentum 
by emitting a hard gluon which can result in the photon energy lying much below 
mb/2. This gluon emission produces a long tail in the photon spectrum below mb/2. 
If E is close to mb/2 a subtle interplay between the two mechanisms takes place. 
As a matter of fact, the very definition of what can be called 'soft gluon' from 
the cloud (the corresponding effect will be described by the primordial distribution 
function) and what can be called 'hard gluon' (the effects due to the hard gluon 
emissions are to be incorporated additionally) is a matter of convention. We must 
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introduce a normalization point /i, and everything softer than will be referred to 
the cloud while everything harder than /z will be treated in perturbation theory. 
Correspondingly, all theoretical quantities - the primordial distribution F{x), A, 
and so on - become /i-dependent. In principle, the final prediction for the physical 
quantities (e.g., the photon spectrum) should be //-independent. Since different 
approximations are made in the soft and hard domains, only an approximate /i 
independence will hold for a physically reasonable choice of /i (see below). 

Technically, the proper way to introduce the primordial distribution function, 
F{x), with 



x = = i^E-^j ' ^ = MB-m,, (5) 

is through consideration of the transition operator. The relevant operator prod- 
uct expansion for the transition operator runs over twists, not dimensions. An 
exhaustive discussion can be found in Refs. |T0[ - [Q; here we remind only that 



theoretically one predicts the moments of F{x) as the B meson expectation values 
of some local operators, see below. The observed photon spectrum is then obtained 
as a convolution of the primordial distribution with the spectrum appearing at the 
level of the perturbative gluon emissions. 



dE y ^ J » 

Integration over y runs from — oo to 1 (more exactly, the lower limit of integration 
is yo = —rriQ/A but we will consistently ignore this difference). One should keep in 
mind that dV^^^*/dE is nonvanishing only in the interval (0,mb/2). We also note 
that one can literally use the perturbative spectrum in the above equation only as 
long as one does not apply it to the very low energy part, ~ A; a more accurate 
expression for this case is discussed in Sect. 5. Both theoretical components in 
Eq. (^ will be considered below in detail. 

Prior to submerging into the discussion of the nonperturbative and perturbative 
distributions, a comment is in order concerning the locality of the weak vertex in 
Eq. (|^). For large mt (we now know that mt lies in the ballpark of 170 GeV) the 
main contribution to the weak vertex comes from distances of order m^^ which 
are indeed short in the scale m^^, so that this part of the vertex is truly local, 
and our consideration is applicable in full. However, some small fraction of the 



amplitude (according to estimates of Ref. [^, less than 10% of the total width) 
comes from the part of the amplitude associated with the loop momenta ^ 
k'^ ^ fTti in the penguin integral. (The corresponding photon is radiated by the low 
momentum virtual c quark in the penguin graph.) Numerically this effect is, thus, 
rather insignificant; but even this insignificant effect can be properly treated within 
the OPE-based approach as long as one is interested in the endpoint spectrum. 

It turns out that the nonlocal part of the decay amplitude can modify the ab- 
solute decay rate (the endpoint domain and the low E part of the spectrum are 
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modified differently, generally speaking); however, the shape of the endpoint spec- 
trum is not affected. This property is the reflection of a universal nature of the 
corrections to the endpoint spectrum which do not depend on the underlying struc- 
ture of the decay amplitude, a feature well-known in the theory of the Sudakov-like 
perturbative corrections. The universality of the nonperturbative effects has been 
discussed in detail in Ref. [0, and the above fact can be demonstrated by following 



a similar line of reasoning. Let us briefly illustrate our point. 

Nonlocality arises due to the gluon emission from 'inside' the vertex (structural 
gluons). In particular, if the internal quarks (in our case the charmed quarks) are 
relatively light, the amplitude at the level of 0{as) corrections becomes complex due 
to real intermediate states that become kinematically possible in the decay. And 
even in this case the relation similar to Eq. (|^) holds. It rests on two key facts: first, 
due to the gauge invariance, the leading twist effects in the transition amplitude 
depend only on the total momentum k of the hadronic system recoiling the photon. 
In the endpoint region k"^ <^ ml and, therefore, the momentum k practically does 
not differ from what one has at the parton level (i.e. transition of the free b quark 
into the free s quark where fc^ = ); in particular, the same light-cone distribution 
function F emerges. 

The second crucial property is that in the kinematics at hand the singularities 
in k"^ associated with the penguin amplitude itself lie far enough from the physical 
point k"^ = mb{mb — 2E) where the transition operator is considered. This ensures 
the following: in considering the transition amplitude in the complex plane at small 
fc^ and performing the operator product expansion (with the aim of obtaining the 
discontinuity of the transition operator) in the way it was done in Ref. |jl2[, one 
picks up only the contribution of the physical cut corresponding to the decay chan- 
nel of interest, 6 — > 7 -|- strange jet . The contributions of other cuts, although 
yielding a complex decay amplitude, are analytical in this kinematical region and 
have no discontinuity. In other words, the presence or absence of structural gluons is 
not essential for the spectrum; only the singularities of the perturbative transition 
amplitude matter here. A more detailed presentation of these arguments will be 
given elsewhere. Here we just note that one gets the same physical spectrum as 
shown in Eq. (^ even for a nonlocal penguin amplitude - the only difference is that 
dTl^^^/dE must be understood now as the perturbative spectrum generated by the 
total penguin amplitude, including its nonlocal part. 

Reiterating, some corrections due to nonlocality of the bsj vertex appear in the 
total width, in the absolute normalization of the endpoint rate and in the shape of 
the lower part of the spectrum. The shape of the endpoint spectrum is given by a 
universal formula. The lower part of the spectrum or the total width, however, do 
not require a particular care - they are insensitive to the primordial distribution, 
nor one needs to carry out summation of the high-order perturbative corrections. 
Therefore, taking account of the nonlocality of the weak vertex here can be done in 
a straightforward manner for both perturbative [|26| and nonperturbative correc- 
tions. 
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3 Ansatz for the Primordial Distribution 



In this section we ignore perturbative gluon emissions and discuss only the soft 
heavy quark distribution function which, in this approximation, is independent on 
the normahzation point. We suggest the following ansatz 

= ^0(1- x)e'=^(l-x)" [1 + 6(1 , -oo<a;<l (7) 

where 

2 



x = = i^E- ^ j , A = Ms - mfe , 

E is the photon energy; the parameters a ,b ,c are positive while k is an integer, 
k = 1,2, ... . We do not consider the integer as a fit parameter; rather several 
independent ansdtze corresponding to different values of k are introduced. We will 
focuss here on the case k = 1, and only briefly comment on modifications for k > 1. 

With this choice of the parameters the function F{x) is positive everywhere in 
the physical domain and has exponential fall-off at large negative x. The latter 
property ensures the existence of all moments [^. Here and in what follows we 



will consistently work in the leading non-trivial approximation in 1/m;,; all terms 
which are additionally suppresed by 1/ m^, are neglected. Then in the absence of the 
perturbative corrections 



TdE A V A 

The primordial distribution function is to be properly normalized. In particular, for 
the first two moments one has 

ao = J dxF{x) = 1 

ai = J dxxF{x) = . (9) 

The first equation determines the normalization factor whereas the second one 
determines the choice of the photon energy reference point to be equal to mb/2. The 
latter can always be satisfied by means of an appropriate substitution 

(1 — x) r(l — x) 

which amounts to the change 

b T% 



TC 



Therefore, our function in fact has two free parameters: the value of a and the ratio 

/3 = b/cK 

In general, our ansatz for F{x) behaves as an arbitrary (non- integer) power of 
1 — X near the kinematical endpoint. This seems to be reasonable since purely 
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perturbative bremsstrahlung effects lead, generally speaking, to a similar behavior, 
as discussed in the next section, and the change in the normalization point leads to 
a change in this power. 

The two next moments of the primordial distribution are given |jl2|] (see also [Q) 

by 

2 

a2 = J dxx^F{x) = , 
a, = l dxx'F{x) = -^^^^{B\Y^ grht%qe^oq\B) (10) 

where in the second equation the sum runs over all light quarks and are SU (3)c 
generators. The matrix element of the kinetic operator is defined as 

Not much is known at present about even these lowest moments of F. The best 
theoretical estimate of /x^ existing today is derived from the QCD sum rules ||30 
and suggests the value 

III ^ 0.5 GeV^ 



see also for an updated discussion). This estimate is in good qualitative agree- 



ment with the model- independent bounds |32, 33, 34] and a semiphenomenological 



estimate of Ref. |5[. For the value of A ^ 500 MeV |0|, |31], |3|] one would ex- 



pect then 02 ~ 0.7. The estimate [T^ of the expectation value of the four-fermion 
operator in the third moment leads to —03 ~ 0.3 which, however, is even more 
uncertain. 

One can readily see that for the function of a reasonable shape it is not easy to 
obtain a large value of the second moment. For example, the maximal value of 02 
for a symmetric positive function which vanishes for |x| > 1 is 1, and if additionally 
it exhibits only one peak then even 02 < 2/3. The actual distribution function does 
not have to (and, in fact, cannot) be symmetric. Nevertheless, it is clear that larger 
values of 02 > 1 can be reached only at a price of having a relatively pronounced 
long tail towards the negative values of x, which, in turn, would imply large negative 
values of the third moment 03 |T^. In fact, strict inequalities can be easily obtained 
for any positive function F 

03 < - - 02 - - ; (11) 




4 V 2 



the first inequality is somewhat stronger than the one obtained in Ref. |]T2 . 

To derive these inequalities one merely observes that for any t the integral over x 
over the function {l — x){x — t)'^F{x) running from —00 to 1 is positive; on the other 
hand this integral is a second order polynomial in t and, therefore, its determinant is 
negative. In a similar manner one can get inequalities incorporating higher moments. 
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In view of poor information available about the moments and a crude nature of 
experimental data on the photon spectrum we mainly limit ourselves to the simplest 
ansatz corresponding to k = 1. The resulting function for k = 1 always has a simple 
one-peak shape. It turns out that within this k = 1 ansatz the value of as is 
strongly correlated with 02, and at fixed 02 a rather limited range of values for 03 
can be reached. For example, with 02 = 0.5 we get 0.47 < —a^ < 0.5. In Fig. 1 
we show the domain of values of a2 and 03 which are accessible within the ansatz 
of Eq. (0) with k = 1 and k = 3. As one can see from the figure, again at the 
same point 02 = 0.5 but for k = 3 the flexibility of the ansatz is much better; in 
fact, at this point 0.29 < —a^ < 0.55 if = 3. The limitations imposed by the 
particular choice of the functional form for the primordial distribution may, under 
certain circumstances, turn out to be important for the phenomenological extraction 
of the hadronic parameters. In particular, if the values of A and /i^ are such that 
they lead to the parameters that lie in the vicinity of the boundary for 02 shown 
in Fig. 1, one may suspect a significant bias imposed by the limited flexibility of 
the model. In this case it will be advantageous to go beyond this ansatz to extend 
the flexibility. This also may be desirable when one has more precise experimental 
data. The simplest way to modify it is by considering larger values of k, say k = 3. 
Then one can easily have less trivially looking primordial distribution functions; in 
particular, it can have two maxima. Correspondingly, the domain of variation of 03 
at fixed 02 becomes much wider which is illustrated by the dotted line in Fig. 1. On 
the other hand, if the result of the experimental fit happens to fall well inside the 
allowed region, one can assume that the bias caused by the concrete choice of the 
distribution function is minimal. 

We pause here to make the following remark. As it will be discussed in detail 
in Sect. 4, the consistent treatment of nonperturbative effects in l/mq expansion 
by means of OPE requires the introduction of a separation scale fi to ensure that 
the domain of gluon momenta below n gets excluded from the coefficient functions. 
In such a way, in particular, the mass of the heavy quark, rriQ that enters the 
expansion is not the pole mass that cannot be consistently defined at the level of 



nonperturbative effects ||3^, |3^, ^ but, rather, the well-defined running mass ^8 



the parameter A is to be understood as A(/i). Although this difference becomes 
relevant only when radiative corrections are considered (they are ignored so far), 
this fact suggests that the proper numerical value of A used in the above numbers 
must be larger by ~ 100 MeV EO] than the estimate based on the one- loop pole 



mass routinely used in the applications of HQET |^ . Then the expected theoretical 
values of the second moment 02 group around 0.4 which seems to be an easy value 
for our distribution function even at k = 1. 
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4 Perturbative Corrections 



In this Section we forget for awhile about nonperturbative aspects and discuss the 
effect of the perturbative QCD corrections. The tree-level spectrum is given by the 
monochromatic photon line at E = m;,/2, Eq. This parton-like formula is 

strongly modified by the gluon bremsstrahlung. The emission of gluons makes the 
energy distribution continuous in the whole kinematical range < E < mh/2 , still 
strongly peaking near the endpoint. We will write it in the form 

TdE dE ^ ' 

where the formfactor S{E) is constant outside the interval < E < mb/2 , 
S = 1 at E <0 and S = at E > mb/2. 

4.1 Sudakov's Double Logarithms 

Let us first approach the problem of the gluon emission in the problem at hand in 
the classic approximation of pure double logs. The inclusive 6 — 57 transition is an 
ideal example since the decaying b quark is at rest while the s quark produced is 
ultrarelativistic; the corresponding drastic rearrangement of the color field manifests 
itself in producing two logarithms for each power of as in the loop integrals. In 
this way we obtain a reasonable qualitative picture of the perturbative spectrum 
by summing up all terms of the form a"ln^" [nib /{nib — 2E)). The result of the 



summation is given by the Sudakov exponent 

5di(E)=e-3^''^'™''/(™''-2^) . (13) 
Correspondingly, in this approximation the photon spectrum reduces to 
Idldi 8a, In (mfe/(mb - 2E)) 

r-JiT = ^ TTE^ bA\[,E) . (14j 

i dE Jvr nib — 2E 

The crucial practical question is what value of as is to be used in this equation. 

In the classical Sudakov approximation the running of a, is not seen since it 
affects only terms beyond the pure double log approximation. Moreover, is also 
not introduced so far. Below we will try to improve this approximation by including, 
at a certain level, the effects of running. For the time being one can use, for example, 
as normalized at the scale mb- It is clear that the double logarithmic spectrum in 
Eq. (|l3|) is quite close to the delta function; the width of the peak is of order of 
nib exp (—const / ^/al) . 

The above double log expressions have rather transparent interpretation which 
allows us to easily improve the precision later on. To this end let us recall that the 
first order probability of emission of the massless gluon in the 6 — >■ 57 decay is 

Vdojd^^ 37rcj(l -cos?9) ^ ' 



10 



where a; is the gluon momentum and is its angle relative to the momentum of the 
s quark. In this expression it is assumed that u <C rrif,. The photon energy in the 
presence of a gluon in the final state is given by 

^ _ ml- 2uj{mb- u) {1- cos 9) ^mt, fcj 

2mfe - 2u;(l - cos^) ~ 2 Au ' ^ ^ 

k± ^ iu-t} ; 

combination of Eq. (^) and Eq. (|16|) reproduces the first order continuum part of 
the photon spectrum obtained by expanding Eq. (JT^) in a^. 

The full Sudakov expression for the formfactor S is obtained in the following 
simple way. One starts with computing the (first order) probability w{E) for the 
gluon to be emitted with such momentum that the photon gets energy below E. 
This probability is obtained by integrating the distribution (p3D with this constraint 
using the kinematical relation (|TB|): 

^■(E) - f ^^^^ 14^, . (17) 



Tdud^^ \Auj \2 J J 371 ruf, - 2E 

The function w{E) has the meaning of the probability of emission of a sufficiently 
hard gluon lowering the photon energy below E. The all-order summation of double 
logs in the Sudakov formfactor amounts then to merely exponentiating this proba- 
bility [0, 

S{E) = e-'"^^) , (18) 

where w{E) is given in Eq. (|1^. As we will see shortly this classical result borrowed 
from QED can be used in QCD only provided E does not approach mb/2 too closely 
4^ . Otherwise, the gluon may become too soft, and such gluons, by definition. 



must be excluded by the perturbative calculations. 



4.2 The Running ag 

The standard Sudakov line of reasoning does not allow one to collect subleading 
terms - those which have higher powers of as for the given number of double log- 
arithms. Needless to say that in QCD it is practically important to determine an 
appropriate scale of as that enters the Sudakov formfactor, which amounts to the 
inclusion of subleading terms. Separation of the subleading terms in the background 
of double logs is a notoriously difficult task. However, one can readily single out a 
class of the subleading terms, namely those which are responsible for the running 
of as, leaving all others aside. A theoretical justification of this procedure is pro- 
vided by the fact that b, the first coefficient in the Gell-Mann-Low function, is a 
numerically large parameter; those graphs which refiect the running of ag do contain 
powers of b while all others, presumably, do not. (In the Coulomb gauge the relevant 
graphs are just the bubble insertions in the gluon line. In arbitrary gauge they can 
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be traced through their specific dependence on the number of the quark flavors). 
The prescription of keeping only this subset of graphs is the heart of the so called 
BLM approach ^5\. In our problem it is natural to combine the BLM prescription 



with the Sudakov-type consideration. 

The original BLM approach was engineered as a scale-setting prescription based 
on the "6-dominance" estimate of the next-to-leading correction. We will first see 
how it works in the problem at hand, and then proceed to a generalized scheme which 
has been used for a long time in theoretical discussions of the renormalon divergences 
of the perturbative series, and was applied recently for numerical improvement of 
the perturbative estimates in QCD |^ ^ . This generalized approach, as it was 
used previously, is applicable to the situations when one deals with a single gluon 
line (dressed by all bubbles). It amounts to inserting the unexpanded expression 
for the running coupling constant ^^(/c^) ~ 47r/(61n (fc^/Aqcj-,)) inside the integrand 
of a one-loop Feynman graph |^ which depends on the gluon momentum k, with 



the subsequent integration over k. In the problem of the Sudakov formfactor one 
sums over arbitrary number of gluon lines. However, since the resulting expression 
is just the exponent of the one gluon emission (i.e. in the given approximation all 
emissions are independent from each other), the application of the generalized BLM 
prescription is quite straightforward and legitimate as long as we stay within the 
approximation of the exponentiation. 

The anatomy of the leading double logarithms presented above allows one to 
easily carry out the corresponding improvement. Indeed, for the emission of a soft 
collinear gluon {uj <^ and -C 1) the effective running coupling enters at the 
scale k±, i.e. as — > as{k±). This fact is well-known in the literature (see, e.g. [pO[] ) 
and it can be substantiated as follows. To determine the relevant momentum scale 
involved, let us introduce a fictitious gluon mass A in the denominator of the gluon 
propagator (we will need this technical trick later on anyway). We then note that 
introducing the gluon mass leads to exactly the same (main) kinematical impact as 
the transverse momentum of the gluon; for instance, the off-shellness of the quark 
propagator is given by k"^ ^ {k\ + \^){mi,/2uj) if the on-shell gluon is emitted in 
the kinematics giving rise to leading logarithms. In other words, the cutoff over the 
spacelike gluon momentum in the rest frame looks like a cutoff over the transverse 
momentum if seen in the infinite momentum frame. 

If so, taking account of the running of one has 



dk\ , k 




frrib 




\ 2 




dk± , 


rrib 


k± 


ki 



In^asik^) + r ^In^asikM (19) 

e •IJein^ k± k± I 



where 

e = rrib — 2E 
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The formfactor S{E) is still given by exp {—w{E)). If one neglects the running of 
then Eq. (|19]) obviously reproduces the standard Sudakov exponent. Eq. (|T^) has 
been displayed in the literature in various forms (see, e.g., [^). 

One can further use the explicit expression ([T9|) to estimate the effective value 
of the strong coupling which enters the Sudakov formfactor. Expanding 

as{k±) = asia) - ^al{a) In — 
Ztx a 

one determines the scale a for which the average of \n{k_<^/a) in the logarithimic 
integral in Eq. (|TU[) vanishes. It is not dificult to check that 

= rrihe . 



As a result we arrive at the following expression 

S{E)=e ^i-'-<>A (20) 

exactly coinciding with the known next-to-leading formula [ pT]] . 

We hasten to add that Eq. pO|) is valid only outside the endpoint domain, 
i.e. when e is parametrically larger than A. (Of course, e has to be parametrically 
smaller than rrib for the logarithmic approximation to have sense.) The Sudakov 
formfactor in Eq. ( pO|) has a peculiar feature. If e ~ Aqcd the strong coupling that 



enters Eq. (pO[) corresponds to the scale a ~ yAQCD^^fe ^ Aqcd- It seems, naively, 
that we are deep inside the perturbative domain. This, however, is not true. As 
soon as e approaches the scale of Aqcd, the integral (p!7| ) defining the exponent of the 
formfactor includes the region of the small gluon momenta (of order of Aqcd) where 
the perturbative consideration is inapplicable. This regime calls for the explicit 
introduction of the normalization point fi and must be treated separatly. In the 
next section we will do a more careful analysis pertinent to the endpoint domain. 

In view of importance of this fact, let us repeat the statement: although the 
expression for the radiative corrections ( PP] ) involves asiy/errib) and therefore, at 
first sight, allows descending down to the heart of the endpoint region, actually 
it is applicable only when e ^ Aqcd (or e > /i = several units x Aqcd in the 
logarithmic approximation) when the integral in Eq. (|19|) runs over the perturbative 



domain. This fact has been noted already in the first of Refs. |£T|. On the other 
hand, as long as e > Aqcd , this asymptotics is correct (such an extravagant "step- 
like applicability" is the usual consequence of the logarithmic approximation). This 



property of the radiative corrections has been used in Ref. p2 |. 

Let us stress that the breakdown of the perturbation theory when E approaches 
the endpoint by the distance of the typical hadronic energy scale signals that non- 
perturbative effects must appear in the spectrum when e decreases down to ~ Aqcd- 
In other words, as it usually occurs the existence of nonperturbative corrections is 
indicated by the perturbation theory itself. 
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4.3 The Endpoint of the Sudakov Formfactor 



The description of the domain e ~ Aqcd (where nonperturbative effects show up) 
requires a more accurate treatment of purely perturbative corrections as well; in 
particular, introduction of the renormalization point n is necessary, following the 
Wilson procedure where all momenta below are excluded from the perturbative 
part of corrections (for a dedicated discussion see Refs. 



The excluded soft 

contributions are referred to the primordial distribution function which, then, also 
explicitly depends on /i; in particular, the local operators that determine its moments 
are normalized at the scale /i [0. It is advantageous to choose the normalization 
scale /i as low as possible still ensuring the perturbative regime as(/i)/7r -C 1. Prac- 
tically, we keep in mind that fi = several units xAqcd- This necessarily modifies 
the standard expressions for the Sudakov formfactor. 

Consistent introduction of the separation of large and small distances to all orders 
in perturbation theory is an unsolved technical problem. In the approximation of 
interest, however, the task can be easily accomplished in different ways. First, one 
may ascribe a fictitious mass in the gluon propagator. This route is very convenient 
even in a purely technical aspect |^ and will be discussed in Sect. 4.5. Here for 
illustrative purposes we will adopt even a simpler method cutting off all integrations 
over the gluon transverse momentum at k± = /i. As a matter of fact, such cut off 
yields the same result for double logarithmic integrals. In different contexts similar 
cut off procedure was used for the same purpose - discarding the soft contribution 
- e.g. in recent work [ p2| . 

Then we arrive at 



w{E;fi) = — 
3tt 



- — In — as{k±) 9{k± - /i) 
k I e 



r^b dk± rUb 
+ I -J — \n — as{k±)9{k±-ij) 



(21) 



where w{E, fi) is the same probability as in Sect. 4.1 with the gluon momentum cut 
off. 

For the one-loop as{k±) the integrals are readily calculated analytically. One 

has 



In 



a 



as[mb) asirub) 
for e > /i ; 



In 



as{e) as{y/erfu) 



(22) 



moreover, 



8 /27r^2 



In 



a.,{y/emb) 



1 a,(/i) 



3tt \ b J [as{mb) as{mb) ^^(/x) as{y/em^) 



Aln/^fl-ln^ii/^ 



(23) 
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tor — < e < /i 

rUb 



and, finally, 



8 /27r\' 



37r V & / 



1 as(/i) 



b , mb 
— In — 

ZTT /i 



(24) 



The photon spectrum is given by 

1 dTP^'^ 



r dE 



for e < — . 



dE 



9{2mb - E)e 



-w(E]ii) 



(25) 



The spectrum implied by Eqs. (|^) - (|25|) has a peculiar feature: w(i?) does not 
depend on e for e < /x^/mb and is equal to w{^] jj) given by Eq. (plj). Therefore, 
differentiation in Eq. (ESf) leaves us with the vanishing spectrum in the interval 



rrib_ 
2 



2mh 



<E< 



2 



plus the delta function corresponding to the two-body decay at the perturbative 
endpoint E = mi,/ 2; the height of the two-body peak is suppressed by the factor 



exp{- 



-w 



mb- 

■ 2 ' 



This behavior of the spectrum is not surprising. Indeed, as was mentioned above, 
introducing a cut off in k± is kinematically equivalent, in our logarithmic approx- 
imation, to ascribing a finite mass n to the gluon. It is quite clear then that the 
gluon mass /i purely kinematically does not allow to fill in the fake window (or mini- 
window) E > Eq^h) = {mb/2) — {2mb)~^n^ the perturbative spectrum. Moreover, 
simultaneously it lifts the absolute suppression of the two-body mode, characteris- 
tic to the classical Sudakov expression. This miniwindow, somewhat aesthetically 
unapplealing feature of the perturbative spectrum obtained above, is not physically 
significant, though, and could be actually eliminated, as discussed below. 

A feature which is, perhaps, more surprising at first sight, is the /i dependence of 
the spectrum in the region /i > e ^ /i^/ (2mb) corresponding to the invariant mass 
squared of the recoiling hadronic system much larger than /i^ . This is a consequence 
of the logarithmically enhanced perturbative corrections involving integration over 
the gluon energies ~ duj/uj. For small uj -C mf, giving the main contribution, the 
kinematical impact of the gluon mass is enhanced: it is given by /i^/2a; rather than 

Below {nib ^ /^)/2 the perturbative photon spectrum does not depend on the 
normalization point, as it should be. 

The theoretical prediction for the perturbative photon spectrum which we present 
here corresponds to consideration of the leading twist operators only. Fine details 
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of the spectrum, at the energy scale smaller than, or of order Aq^j-i/mf,, are shaped 
by high twist operators completely ignored in our analysis. As a matter of fact, in 
order to resolve the fine details of the spectrum at this scale one needs to analyse an 
infinite chain of corrections of all possible twists 0], a task obviously going beyond 
possibilities of the present-day theory. For reasons explained above we push /x as low 
as possible, in practice, a few units xAqcd- Then, at our level of accuracy we can 
not distinguish between the spectrum presented above and the one where the elastic 
spike is moved to the left so as to close miniwindow, or just smeared evenly over the 
miniwindow. The corresponding changes in the moments / x"'F{x)dx are of order 
l/rrij,. As long as one keeps the integrals intact to order 0(m°) one can distort w{e) 
compared to the formula (plf) at will. In particular, one of the ways of eliminating 
the elastic spike at e = is as follows. At e <^ Aqcd one can use arbitrary (growing) 
w{e) and not necessarily the one given by Eqs. (p3D, (|^); assigning w{e) — > oo at 
e ^ smears the elastic peak completely. The prediction for the physical photon 
spectrum remains intact at the level of accuracy we work at in this paper. 



4.4 The Endpoint Domain and the b Quark Mass 

Specifying the theory above, we normalized all quantities including the weak 657 
coupling h, as and at some scale. For the inclusive b decays the appropriate 
scale is ~ m?,. Thus, we start from a high scale mass, say, mb{mb). On the other 
hand, practical perturbative calculations are done similar to QED analysis which 
is phrased, of course, in terms of the on-shell mass. Literal application of our 
expressions, therefore, implies the use of the pole mass, in QCD. This mass 

is defined to any finite order in perturbation theory; in fact, it is a meaningful 
theoretical notion in the theory where a strong coupling regime does not occur. 
Then it is tempting to say that the endpoint of the spectrum is shaped by this pole 
mass, which is higher than m;,(mf,). 



However, it is well known [p7| , pq , 3£] that the true pole mass cannot be consis- 
tently defined in QCD with the precision we are interested in, ^ Aqcd- This exactly 
corresponds to the circumstance explained above in detail, that one cannot compute 
the perturbative spectrum for e ~ Aqcd naively, without introducing the normahza- 
tion point /i. Of course, the problem does not show up in the expressions obtained 
in the pure double logarithmic approximation where the running of a<j is discarded, 
a "frozen" value of as is used and fi is not introduced. In this approximation there 
is no difference with QED, and the corresponding value of rrib essentially coincides 
with the "one- loop pole mass" discussed in Refs. ||3^, |3^, Q . 

We go beyond this approximation, however; correspondingly, the question arises 
as to what mass enters Eqs. (|T^) - (|T^) and below. It is clear that the mass we 
deal with in the endpoint domain in our approach is the would be pole mass in the 
theory where the gluon exchanges with \k\ < are excluded, i.e. the running mass 
normalized at /i, mb{fi). It is important to notice that mb{fi) has nothing to do 
with the dimensionally regularized mass mql/j,) (the latter mass is a well defined 
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theoretical notion but it is irrelevant for OPE in the infrared region /i <C fTiQ). The 
appearence of the running mass normalized at n in all expressions for the spectrum in 
the endpoint domain can be exactly proven in the BLM approximation where there 
are no technical problems with defining the explicit renormalization prescription. 
We have tacitly assumed the BLM approximation in the calculations above. For 
a related discussion see Ref. ||5^. The occurence of justifies the concluding 

remark of Sect. 3. 

The perturbative spectrum does not extend above E = and the window 

between mb/2 and Mb/2 is filled when nonperturbative corrections are accounted 
for. However, in the OPE-based approach the b quark mass, and, therefore, the 
endpoint of the perturbative spectrum, depends on /i. Let us dwell on this aspect 
in more detail. The dependence of the heavy quark mass on /i at /i -C mq is given 
by il,!! 

djj, IT 6 

and, therefore, mq increases with fi decreasing. Calculating perturbative corrections 
we must stop at some /iq ^ Aqcd- Then the mass we "see" in the endpoint region is 
mb{fio) and, in particular, the perturbative spectrum spans up to -Emax = ?7^6(/^o)/2 
(we assume that /xq ~ a few units x Aqcd and consistently neglect the corrections 
~ /ig/2mfe). Assume we can go perturbatively a little bit further, down to a smaller 
value of /i. The perturbative physics below fiQ will populate the spectrum up to 
mfe(/i)/2 and, therefore we will observe the decay events above the original pertur- 
bative "quark endpoint" mb(/io)/2. 

Formally, the whole physics below /xq is now treated in the fio/^b expansion 
and is accounted for by a 'soft' distribution function normalized at the scale /iq. In 
particular, it refers now as well to the region between /i and /io which is still amenable 
to the perturbative treatment. Moreover, one can even introduce a perturbative 
analog of A, 



A = mb{fl) - nibiflo) = Crri rfyU ~ Cm {no-fi) + ... (27) 

Jfj. 71 TT 



Clearly one cannot, though, push /i too low to get an "estimate" of the whole window 

A(/io) = Mb- rribifxo)- 

This remark shows once again that many qualitative features of the nonper- 
turbative dynamics have their seeds already in the perturbative expansion itself. 
Moreover, in QCD it is impossible to distinguish the perturbative effects from the 
nonperturbative ones, and one should rather speak about separation of the short 
distance effects from the long distance ones. 

Eqs. (|22|), (P^l) determine the asymptotic suppression of the domain mb ^ e ^ 
Aqcd when the heavy quark mass goes to infinity, 

■§r In "■'^"^ +const const(e) — Inln 



S^mn ' . (28) 
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It, thus, formally decreases faster than any finite power of tjiq. More instructive is 
the behavior of the spectrum as a function of e when uiq is asymptotically large but 
fixed. Eq. ( PB| ) shows that the formfactor S at small e < /i behaves like some power 

of 



The critical exponent eo actually depends on the normalization point, rather weakly 
though, as long as one stays within the domain of validity of the perturbative treat- 



ment [54]. The variation of eo upon changing normalization point /i is calculable 



and is given by 

^de^ = + Oialif^)) + OM^)) . (30) 

Again, one cannot use this expression with too low of a normalization point, /i ~ 
^QCD, and, thus, determine a universal scaling law for the perturbative distribution 
function. 

The results which were presented in Sect. 4.3 take into account all terms of 

the form (^a^ln^mb/e^ ■ (6a<j/7r)'. In other words we allow to loosen the Sudakov 
logarithms provided that the correction contains b, the first coefficient in the Gell- 
Mann-Low function. Thus, the approximation we work with can be called all order 
BLM-improved double logarithmic approximation (it will be referred to below as 
'natural'). It is similar to the approximations, often rather sophisticated, used in 
applications to problems where the Sudakov-type effects appear. The most crucial 
difference is the explicit introduction of the separation scale /x. Although b is not a 
true free parameter in QCD, its numerical value is rather large, and similar BLM- 
type approximations often work rather well. What is even more important for us, 
within the natural approximation we are able to capture main features of the full 
result. Some of them are of a rather subtle nature and are often misunderstood in 
the literature. 

The BLM prescription treating 6 as a formal parameter allows one to determine 
the appropriate scheme which must be used for the strong coupling - Eqs. (0), 



(pl|), (|22| ) - (^) and others are valid if as is understood in the V scheme. On 
the other hand, say, the precise kinematical boundaries of integrals in terms of e, 
H and rrib cannot be fixed in this way. For example, changing everywhere in the 
integrals determining w{E) e —>■ 2e or — > 2m;, does not lead to the terms of the 

form (^aglTi^ TTLb/ej ■ (foas/vr)' summed up within the natural approximation; their 
consistent computation requires an honest calculation of the next-to-leading terms 
and is very involved. 

An improvement of the Sudakov-type calculations was considered in the recent 



publication |Q. The results obtained there provide one with the better accuracy 
in determining the dependence of the spectrum for asymptotically large mass. 
They can not be directly applied, however, to the solution of our problem, since 
the separation scale fi was not explicitly introduced . The moments |^ we would 
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need to consider require integrating the anomalous dimensions in the region below 
the infrared pole. On the other hand, introducing the normalization point n the 
way we do would violate the scale invariance, a crucial component of the method of 



Ref. 55 



Even leaving this aspect aside, determining the spectrum from the moments 



presented in |^5[ is a difficult problem. In the next section we propose an alternative 
approach which can be carried through rather straightforwardly and is expected to 
provide sufficiently good accuracy in 6 ^ S7 transitions. 

4.5 A Suggestion for Further Improving dVP^'^^/dE 

Combining the BLM prescription with the classical Sudakov calculation is already 
an improvement, beyond any doubt. However, as we have just noted, this natural 
approximation does not distinguish, say, between rrib and 2mb in the argument of 
the logarithms. The logarithm is far from being asymptotically large in the b decays, 
and significant progress in experimental data on the — > s + 7 spectrum will call 
for a more refined theoretical description of the perturbative part of the spectrum. 
One will have to adress two different questions, namely further improvement of 
the bremmstrahlung spectrum, and calculation of corrections to the weak decay 
vertex itself (which is not exactly local). The latter complication has been briefiy 
addressed in Sect. 2, and we focus now on further improving the accuracy of the 
bremsstrahlung spectrum. The consistent summation of all subleading logarithms is 
clearly an unrealistic task, the more so that we want to do it in a way compliant with 
the Wilson OPE, i.e. discarding the infrared domain in the perturbative calculation. 
Therefore, we are putting forward a suggestion resulting in what we will call "the 
advanced perturbative spectrum" (APS), combining the exact first-order result, the 
Sudakov summation of the logarithms and the extended BLM improvement. 

(i) Since In (mfe//i) is not a particularly large parameter we believe that it is 
sufficient to exactly calculate a few first orders in (in the simplest version - 
merely the exact one-gluon expression) still continuing the summation of exploding 
double logarithms. For rub — 4.8 GeV the terms 0{as) seem to be more important 
numerically than even the next-to-leading logarithm, O {{a'^/ir) In {rrib / fi)) (at least 
when the latter is not related to the running of the coupling). 

(ii) One should continue to use the expression for the perturbative spectrum 

(0), 

with the exponentiated form of the Sudakov formfactor (p!8|); however, the exponent 
w{E; ix) is calculated now as the exact first order probability of the gluon emission 
within the kinematical constraint that the resulting photon energy is less than E, 
rather than the double logarithmic probability as in Eq. (p!?!). 

The parameter /i must be introduced in such a way as to isolate and suppress the 
infrared part. In the double logarithmic approximation this goal was easily achieved 
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by discarding the low k± part of the k integration. Now we can do the same job by 
introducing a fictitious "gluon mass" . The density matrix in the gluon propagator 
is kept the same (e.g., in the Feynman gauge it is 60/3) while the Green function 
(A;^ — This "gluon mass" does not lead to any bad consequences in the 
one-gluon graphs (exponentiation does not count), which are the same in the abelian 
and non-abelian theories; in the abelian theory the mass term does not violate the 
conservation of the (abelian) current. The same trick was proposed recently in a 
different context ^ and we refer the reader to these publications for further 

comments. 

To the first order in a, the probability of the gluon emission (leading to the 
photon energy less than E) is given by 



w{E; /i) = asW{E; /j,^ 

11 .a{i-2e:)-i 

e{Eo-E) dE' dps\A{E',p„^)\' (32) 

JO J A 



where 



Eo = — -- — , cos?9=l + — 1 , A = Eq-E. 

2 2mb E' \ ps J 

In the above equation E' is the photon energy, ps is the energy (momentum) of 
the strange quark and i) is the angle between the momenta of the photon and the 
quark; \A{E',ps,^)\'^ denotes the amplitude squared summed over polarizations of 
all particles in the final state and color configurations; it depends on /i through 
kinematics. It is directly calculable but the analytic expression is too lengthy and is 
not important for our discussion. We have introduced above W{E,fi), the probabil- 
ity of having the photon with energy less than E, in the one-gluon approximation, 
with the factor as separated out. This quantity will be useful in what follows. 

The total width r(/i) is obtained by the exact calculation of the decay width 
within the same {viz., first order in a^) approximation The formula for the 



perturbative spectrum obtained in this way has obvious advantages. On the one 
hand, it gives the exact expression for the one-loop spectrum (with massive gluon) if 
expanded in a.^- On the other, it automatically leads to the summation of the double 
logarithms near the endpoint region where the one-loop expressions explode for large 
rrih/fi. At the same time, it respects the exact relation between the contributions of 
real and virtual gluons. This completes the first part of our suggestion. 

(iii) The next step of our program that naturally follows from analysis in Sect. 4.3 
is to BLM-improve this exact one-loop formula by using the running strong coupling, 
thus promoting it to the status of the advanced perturbative spectrum. Again, both 
the exponent w{E,^) and the overall normalization r(/i) must be calculated within 
the same procedure. 

Warning: we do not claim that the exponentiation of the exact 0{as) amplitude 
(or even of the BLM-improved one) reproduces the full result in all orders, as is 
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the case with the standard Sudakov double logarithms. Moreover, we are unaware 
of any adjustable parameters that will control the accuracy of this prescription. 
Nevertheless, arguments can be given that numerically it must have a good accuracy 
in the problem at hand. It is important that the corrections to exponentiation are 
proportional to higher powers of as(/i) (or even evaluated at higher scales), and, 
therefore, are to be calculated completely within perturbation theory. 

Calculation of w{E] /i) in the BLM approximation accounting for the running of 
as can be easily done explicitly; it can be written as the integral over the fictitious 
gluon mass of the width introduced in Eq. (|32D , simple three-dimensional 

integral. The concrete expression is too cumbersome due to lengthy expression for 
the tree-level decay amplitude itself, and will be given elsewhere. The main technical 
complication is that this BLM-improved calculation must be done respecting the 
separation of low and high momentum scales required by OPE. The purpose of the 
discussion below is to show explicitly how it can be accomplished in a straightforward 
way. 

4.6 Outlining a Method for Introducing fi beyond the Lead- 
ing Logarithmic Approximation 

The extended BLM-improvement of the one-gluon graphs per se is a rather simple 
technical exercise for Euclidean quantities where one merely uses the gluon propaga- 
tor in the form as(/c^) / fc^ with the running inside the integrand |5^, ^ 



Moreover, it can be directly generalized for the quantities of the type of decay widths 
which are formulated in Minkowski space ^ ; the results, however, cannot 



be used in the all-order resummed form in the both cases because one encounters the 
unphysical infrared pole in the gluon propagator which has wrong analytical prop- 
erties. This problem is absent in the OPE-based approach where the contribution 
of low momentum domain is excluded from the perturbative coefficient functions. 
Introduction of the hard infrared cutoff is rather trivial in the extended BLM 



approximation for the Euclidean integrals (see Refs. [38, 48, 53 1 for examples). How- 
ever, such method is not directly applicable to the Minkowski integrals. Moreover, 
application of OPE to the inclusive widths introduces some peculiarity absent in the 



analysis of Ref . , and we briefly discuss it here. Details of the derivation and the 



complete discussion of the approach will be given in a separate publication [pOj] . 

Thus, our task is to eliminate the soft gluon part in the coefficient functions. 
This task can be achieved in two steps. First, we soften the theory in the infrared 
domain by introducing the gluon mass. In the BLM approximation the mass term 
can be added by hand; otherwise one could have introduced the gluon mass through 
the Higgs mechanism. If this mass is larger than Aqcd the effective color interaction 
never becomes strong. Then we calculate the coefficient functions in this softened 
theory. 

Even though the theory is now softened in the infrared domain, our result for the 
coefficient function will still contain some (power suppressed) residual contributions 
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coming from the infrared domain. We still have to subtract them, but now the 
subtraction is much easier to do than in the unmodified QCD because the modified 
theory is not strongly coupled in the infrared domain. The subtraction can be carried 
out to any given order, and the procedure will presumably converge. Implementing 
this strategy which seems to be applicable rather universally (at least in principle), 
we go beyond the practical version of OPE pT| . 

A subtle point should be noted immediately. Modifying the theory in the infrared 
domain for calculational purposes we must ensure that it is not modified at large 
momenta - otherwise, instead of solving the technical problem of building the correct 
OPE in QCD we will be merely dealing with a different theory. In the procedure we 
suggest this latter requirement can be met to any given (fixed) order in the power 
series. 

The advantage of introducing the separation scale /i in the above way over the 
hard cut off in the Euclidean integrals, with the subsequent analytic continuation to 
the Minkowski domain, is rather obvious, even leaving aside purely computational 
aspect. In doing so we do not violate the analytical properties of the theory (and, 
therefore, the perturbative coefficients themselves) - the feature which cannot be 
achieved if eliminating of the low-energy contribution is carried out in the "hard" 
way. It is worth emphasizing that the procedure of softening is a technical device for 
OPE-based calculations. In particular, the short-distance parts of the amplitudes 
are not required to be unitary - the unitarity holds only for the complete physical 
amphtude. 

How is this general strategy is implemented within the extended BLM routine? 
The substitution as/k'^ as{k'^)/k'^ for the renormalon chain in the integrand leads 
to the Landau pole at the (Euclidean) point Aqqj-, where no singularities are allowed. 
However, one can merely shift this pole, by hand, to a Minkowski point — /i^. 

Here we use a convention according to which k'^ is Euclidean The one-loop expression 
for the running is implied. Here and below all couplings and Aqcd refer to the V 
scheme. Of course, the large k"^ asymptotics is changed, but only at the level k~'^, 
which is important for power corrections. This drawback can be eliminated through 
any given order in the power series. 

If we want to keep (a finite number) of the power corrections we must ensure 
that by modifying as{k'^)/k'^ in the infrared domain we do not spoil the ultraviolet 
asymptotics, up to the level of the power terms we are interested in. Then we could 
choose, for instance, the infrared regularization as follows: 



^sik^), , r^BLMn.2..A _ , (as{k^) An 1 ^An_ 

^QCD b + ^ 
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47r ' (/i' + A^, 



(34) 

where / is a finite number corresponding to the number of Pauh-Villars subtrac- 
tions. This function coincides with the unsoftened resummed propagator as{k'^)/k'^ 
at /c^ ^ /i^ up to terms ~ yU^'+^/fc^'^'* and, thus, allows one to treat properly (within 
the accuracy of extended BLM) the condensate corrections through dimension 21 
(excluding the dimension of the heavy quark fields themselves). For instance, to 
address the effect of the kinetic energy operator it is a priori necessary to use / > 1 . 
(^BLM^^2. j^g^g Landau singularity in the Euclidean domain and is regular and 
suppressed at A;^ < yU^; in particular the effective coupling vanishes at fc^ ^ 0. 

If this infrared regularization method is accepted, the result of all-order BLM- 
summation for the generic inclusive width, w^^^{fi), (in our case it is w^^^{E; n)) 
can be expressed as an integral over the "gluon mass" over the one-loop value of w 
obtained with the gluon mass A, WlX"^) with a certain weight function 0^'"'^: 

^^"^(/^) = f ^ W{X') r'^'^iX'; f,) . (35) 

This expression is similar to the one obtained in Refs. however, (p^^^ depends 
explicitly on the normalization point fi and is not a function of only the "effective 
gluon mass" A as was the case in the analysis of Refs. [0. Moreover, it has quite 
different properties. It is given by the imaginary part of the gluon propagator in 
Eq. (0), 

</>"^^(A^;/^) = ^ImGL^^(AV) = 

b In^-r^ h7r2 b b n\ 

The first (regular) term in 0^^^, when substituted in Eq. (^5]), represents the 
perturbative contribution of the continuum existing for all invariant masses A^. The 
sign of this continuum contribution is negative; this is a reflection of the asymptotic 
freedom. The delta function terms have the right (positive) sign. Hence, this con- 
tinuum contribution is a correction diminishing the probability. The 5-like terms 
in Eq. ( P^D have the meaning of the probability of emission of the gluon with the 
mass fi (with small corrections ~ yU^" making up for the short distance dependence 
of this width compared to the actual massless case). The appearance of such local 
terms with higher derivatives is typical when one applies the OPE procedure to the 
spectral densities in Minkowski space |^ . 



In principle, derivation of Eq. ( |35D assumes that ^(A^) is an analytical function 
of A^ at A^ <^ yu^. This is always valid in the consistent application of OPE: a 
closer singularity would imply the presence of other "soft" propagators and would 
mean that they have not been properly regularized in the infrared by adding the 
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corresponding operators in the OPE. In our particular case of the first order gluon 
emission probabihty W{E;X), a singularity in the gluon mass A occurs only at 
2mb — ~ fi'^/iTLf,, i.e. only in the region of the "miniwindow" . Therefore, this 
requirement is met. 

In the formal limit, fi"^ —>■ 0, one essentially reproduces the expressions of Refs. |^ 
where in the latter the principal value prescription is used for treating the Landau 
singularity. This limit can not be taken, however, if one wants to make OPE- 
compatible analysis and to address nonperturbative effects. 

The specific form of the infrared softening is not very important. We have done 
it by adding 5-functions in (f)^^^; alternatively, one could have used some smeared 
functions. The only constraint to be observed is that the first / moments of (f)^^^ 
must remain the same. For example, one can consider 0^lm completely vanishing 
in the vicinity of = which would yield the propagator analytical ai = 0. 
Generically, the gluon Green function in the Euclidean becomes then 



G^j^''ik'■,^^) = 5^, 



oo rfA2 0f^(A2;/i) 



A2 



P + A2 



(37) 



As long as the moments are preserved all expressions of this type are equally suited 
for the OPE. Different choices it merely correspond to somewhat different cutoff 
procedure, and the matrix elements of the dimension n operators then differ by 
perturbatively calculable terms ~ (a;s(yu)/7r)'^yu". 

Similar to the approach of Refs. upon integrating by parts the integral over 
the the gluon mass in Eq. ( pSj ) can be rewritten in the form resembling integration 
over the gluon virtuality with an effective strong coupling af^. The quantity that 
plays the role of the distribution over the gluon virtuality A^ is — dVr(A^)/dA^; 
however, af^ depends on both /i and A and also contains local terms: 



w 



BLM 



(E,/i) 



dA' af (A;/i) 



dA2 



+ ^(m/^Vi^(o)) + ^E 

" " n=l 



where 



«f(A;/i) 



2 A;2 b 





(d/i2)'^ 


71 

arctan 




2 





(38) 



(39) 



with af^ (oo) = and a^^(O) = 47r/6. One can formally introduce the "full" coupling 
af^ which combines all terms: 
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b ^1 n\ 

which is legitimate for smooth Vr(A^); otherwise one can use a regular "smeared" 
representation for local terms as was mentioned above. 

For high momenta, when /i^ ^ ^qcd? ^^e effective couplings a'f^X-jfi) 

and ^^''^(A; /i) coincide with the usual strong coupling as{\). However, a'f vanishes 
for small gluon momenta which reflects the fact of excluding the low momentum 
(Euclidean) region instead of assuming a model behavior of the gluon propagator 
in Ref. |^^. One still has to keep in mind that all effective couplings introduced 
above depend (in the infrared region) on fi and, in particular, on the exact way 
the normalization point is introduced; on the other hand, the physical results are 
scheme-independent, and therefore the effective couplings cannot have a direct phys- 
ical meaning in the low energy domain whatever technical convenience they grant. 
Further details will be available from [ |60[] . 



Let us briefly illustrate how relation (^) works in simple cases. Consider, for 
example, the total inclusive width of a heavy flavor, say, the process t ^ b + 
W^(+ gluon) The one loop correction to the decay width r'^^^(A^) decreases 

fast for A^ ^ mf (it does not vanish for A > mj — because virtual corrections 
are present for arbitrary finite gluon mass). On the other hand, for A^ -C {rrit — 
MwY it is practically independent of A. Therefore with the logarithmic accuracy it 
can be approximated by the step-function, r'^^)(A^) ~ 6{m1 — A^)r'^^)(0); then one 
immediately recovers the obvious result that the proper momentum scale for the 
strong coupling in the one loop correction is of the order of 

To summarize this Section, we believe that the general method suggested above 
can provide one with a very good theoretical accuracy of describing the perturbative 
corrections to spectrum in a number of applications (6 — s + 7, semileptonic decays 
etc.) in the OPE compliant way. Let us repeat the necessary steps. First, one 
calculates the decay rate of interest to the first order in with the nonvanishing 
gluon mass. (In the case of the photon spectrum in 6 — s + 7 one calculates the first 
order probability to have photon with energy below E). Then, fixing the normal- 
ization point fi and using Eqs. (P5D-(p]|) one sums up all perturbative corrections 
in the BLM approximation to get the "BLM-improved" Sudakov exponent and the 
total width which explicitly depends on the normalization point through the effec- 
tive coupling. Then the decay distribution is given by Eq. (|31| ) which exponentiates 
the gluon exchanges. Though looking somewhat cumbersome, this procedure can 
be done numerically in a straightforward manner. 

The very same startegy can be literally applied to semileptonic decays as well; 
the only technical difference is that here one first needs to fix the invariant mass of 
leptons and, at the final stage, integrate over q"^. The contribution of the region of 
maximal in the case ofb^u transitions can be treated within OPE as described 
in detail in Ref. . 
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5 Full 6 ^ s + 7 spectrum 

The full observable spectrum is obtained by convoluting the perturbative spectrum 
with the primordial distribution function responsible for the heavy quark motion 
due to the soft gluonic medium. The precise form of the convolution deserves a 
special discussion. 

An obvious similarity exists between the classic theory of deep inelastic scattering 
(DIS) [^] and the OPE-based theory of the spectra in the inclusive heavy quark 
decays |Tl], |12[ . In both cases the expansion runs over twists, not dimensions, and 
the predictions are originally formulated for the moments of certain functions and 
are given by matrix elemets of local operators times calculable coefficient functions 
(the latter carry all dependence on the large parameter, in DIS or mq in the 
heavy flavor decays). The operators must be normalized at a scale which stays finite 
when or tuq tend to infinity. Then the statement of factorization of the large 
and short distance contributions in the moments is translated into the convolution 
of the primordial soft function with the one describing perturbative evolution due to 
the hard gluon emissions. In DIS this is the essence of the Gribov-Lipatov-Altarelli- 
Parisi equations . 



However, there are some differences between the two cases - in essence, kinemat- 
ical - and a straightforward extention of the parton-model language and all formulae 
we are used to in DIS to the theory of the spectra in the inclusive heavy quark de- 
cays can be unjustified. The main difference is as follows. The moments predicted 
in DIS are the moments of the structure functions -F(x), < x < 1, with respect 
to the Bjorken variable x. If we know a large but finite number of moments (i.e. a 
number not scaling with when oo) we know F{x) in the whole physical 

domain < a: < 1 (with a certain accuracy). Now, in the problem of spectra the 
appropriate variable x is defined in Eq. (^). The knowledge of a finite number of 
moments with respect to this x determines the distribution function only in the 
window < X < 1 and in the nearby endpoint domain — 1 ^ x < 0. One needs to 
consider the moments through order > mg/AQCD to address the whole spectrum, 
including the low E domain, in the heavy fiavor decays. For such high moments the 
standard analysis does not apply. Thus, the standard convolution of the primordial 
and hard components (e.g. Eq. (79) in Ref. 0]) is valid only in the endpoint 
domain; what form is valid for lower values of is a question to be answered. 



The reasoning in Ref. was applied to the endpoint spectrum. In the present 
paper, for practical purposes, we need to have an expression in the whole kinematical 
region < < Mb/2. Fortunately, it is very easy to answer the question what hap- 
pens outside the endpoint domain. We know for a fact, from the OPE analysis, that 
at these values of E the motion of the heavy quark inside B affects the physical spec- 
trum only at the level of corrections. Since such corrections are consistently 
neglected in the present paper the proper convolution formula should be written in 
such a way that outside the endpoint domain the physical and perturbative spectra 
be the same up to terms 0{m^). 
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In fact, the answer can be prompted by a physical picture where one takes the b 
quark at rest and then boosts it to a moving reference frame to take account of the 
primordial motion. This suggests that 



where 



E' = E- -—X 



2 Mb 

This expression is obviously equivalent to Eq. if E is close to nib/ 2, up to higher 
orders in l/mft distinguishing Mb from m^. Moreover, both forms of the convolutions 
are always equivalent where the perturbative spectrum is smooth and can be Taylor 
expanded. Say, outside the endpoint domain (but not for E ~ Aqcd) the difference 
between evaluating the perturbative spectrum at E, E' or E — Ax/2 is noticeable 
only in the next order in l/rUb (and for the convolution, only in 1/ml) which is 
precisely what we wanted. 

The actual difference between Eqs. (|) and (^Tj) can show up only at the points 
where the perturbative spectrum is singular. In principle, it hapens at two points - 
ai E = mb/2 (for the massless s quark) and ai E = 0. At E = mb/2 the formulae 
coincide, as was mentioned above. One must, therefore, address only the point = 
where the perturbative spectrum contains explicit 0{E) just due to the photon phase 
space. 

Before doing it, let us comment on Eq. (l4l| ) from a different perspective. This 
equation sums up all leading twist contributions in the expansion of the transition 
amplitude, which shape the endpoint spectrum (for decays into light quarks) and is 
formally valid everywhere, < -E < Mb/2. However, outside the endpoint domain 
the summed terms are not dominant and give corrections of the same order as some 
of the higher twist operators neglected in the expansion. 

Now, consider small E. A formal derivation of the full spectrum can be carried 
out using the approach of Ref. where the complete transition amplitude is con- 



sidered first in the complex plain. One then observes that originally, instead of the 
perturbative spectrum, the convolution goes with a certain structure function (the 
discontinuity of the hadronic tensor corresponding to the quark process discribing 
the hard gluon emission). This structure function is directly related, up to trivial 
kinematical factors, to the perturbative spectrum in the physical domain, but it is 
smooth at -E = 0; the vanishing of the perturbative spectrum at < is ensured 
by a kinematical factor representing the photon phase space which contains 0{E) . 
This structure function has singularity only near E = mb/2 and, perhaps, at some 
other distant unphysical points. 

Since the structure function above is smooth at zero, the convolution with the 
primordial distribution function in the small E domain has no effect whatsoever; 
one just reproduces the perturbative result. This is not what is done in Eq. (^ 
which merely ignores the phase space factor singularity absent for other values of 
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E. Taking account of this subtlety at small E leads to Eq. (^TJ) which is written in 
such a way that this effect of the primordial smearing is just absent at all for small 
E: the smearing goes over the interval in E which is itself proportional to E; more 
essentially, the interval never covers the point E = 0. The last factor in Eq. is 
introduced ad hoc to ensure the correct normalization (to unity); it is formally of 
higher order in 1/mi,. 



6 Analysis 

At the present level of experimental data it seems premature to use APS in the 
numerical analysis. Thus, we return to the approximation called "natural" in Sect. 
4.4 - exponentiated double logarithms all order BLM-improved. In other words, in 
practice we use the expressions for the perturbative spectrum given by Eq. (^Sf) with 
the Sudakov exponent in Eqs. (p^)-(]2^). 

The form of the convolution exploited to obtain the theoretical spectrum of direct 
photons in 6 ^ s + 7 decay is given in Eq. (|4l| ) . The value of the strong coupling is 
taken following from the one loop value Aq^j-, = 300 MeV (in the V scheme) which 

nearly corresponds to q;^(2.3 GeV) = ^^''^'^(l GeV) = 0.336 as suggested by the most 
accurate low energy analysis ||6^. In view of the poor experimental data available 
we limit ourselves to the primordial distribution functions with k = 1 only. 

Fig. 2 shows the full spectrum (solid line) obtained for the theoretically preferable 
choice of parameters /i^ = 0.5 GeV^ |Q |31| and mfe(/i) = 4.67 GeV 0, ^ where 
we use the normalization scale /i = 0.7 GeV; the dependence of rrih on the normal- 
ization point is assumed to be given by Eq. (p6D |6^ . The dashed line shows the 
corresponding nonperturbative spectrum which thus is based on A(/i) = 610 MeV 
(02 = 0.45 obtained with a = 0.4 and j3 = h/c = 6.4, see Sect. 3). We also show 
purely perturbative spectrum stemming from Eqs. (p^)-(p5D with Aq^d = 300 MeV 
and /i = 0.7 GeV; the area under the ^-function peak constitutes numerically 0.45 of 
the total perturbative width. It is clearly seen that the account for the perturbative 
corrections essentially modifies the resulting spectrum [|12[: first of all, it decreases 
the height of the maximum by a factor of 2.5, and also noticeably broadens it. 

The exact position of the normalization scale /i = 0.7 GeV used above was chosen 
rather arbitrarily, just to meet best the two usual requirements - to have it as small 
as possible compared to and, on the other hand, still have the perturbative ex- 
pansion parameter ^^(/i) / vr sufficiently small. Certainly, this still allows a significant 
variation of this scale. Fig. 3 illustrates the dependence of the full spectrum on the 
choice of /i; however, varying /i we must supplement it by related change in rribdj,) 
in the perturbative spectrum and A(/i) in the primordial distribution function. We 
used the literal form of Eq. (pBD (see also Ref. [0]) to determine the corresponding 
masses. In principle, in the test of the /i independence of the physical spectrum the 
value of /i^ should also have been varied as a function of normalization point. We 
did not do that since numerically the most important effect is due to the change of 



28 



A(/.). 

The solid line in Fig. 3 is obtained with the same, central, choice of parameters 
as in Fig. 2; dashed line represents the high normalization point /i = 1 GeV and 
the long-dashed curve shows the case of low fi = 0.4 GeV; the overall normalization 
is the same for the three curves, unity. If mfe(0.7GeV) = 4.67 GeV then the cor- 
responding mass values for other normalization points are mb(lGeV) = 4.58 GeV 
and mb(0.4GeV) = 4.83 GeV. The value of /i^ = 0.5 GeV^ was kept fixed which 
assumed rather different values of 02, 0.34 and 0.83, respectively. The normalization 
point /i = 0.4 GeV is probably too low to be taken seriously. It is considered only 
for illustrative purposes, as an extreme case. 

Although both the perturbative spectrum and the primordial distribution sep- 
arately look very different for the three choices of the normalization point, the full 
spectra are very much alike. We stress once more that it was crucial to use the 
proper normalization point dependence of the heavy quark mass in order to get this 
similarity; using a universal, or "true", mass as sometimes suggested would lead to 
essentially different spectra. This transparent demonstration of the physical rele- 
vance of the running b quark mass is an exact analog of a similar example analysed 
in Ref. ||3^ for the case of the small velocity (SV) kinematics. 

The three curves in Fig. 3 are not exactly identical. A small difference emerges 
for a few reasons. First, it reflects the finite value of the b quark mass; therefore 
higher terms in l/rrib neglected in our analysis give finite contribution. We also 
used approximate expression for the perturbative spectrum. Second, it is due to 
the rather limited flexibility of the used ansatz for the primordial distribution: in 
principle, the functional form we adopted is not preserved by the renormalization 
from one value of fi to another. This makes it important to use a physically justified 
choice of /i. However, the very fact that the full spectrum does appear to depend 
weakly on fi indicates that our approximations work quite well. 

Turning to comparison of theoretical predictions with data we have to admit 
that the latter are rather crude at the moment to draw any well justified definite 
conclusions. Nevertheless it is fair to say that the theoretical spectrum for the 
preferred values of parameters fits well with what is seen in experiment. In any 
case, keeping in mind the fact that the experimental analysis which lead to the few 
experimental points we used, was rather involved, their naive fit with any theoretical 
prediction per se is hardly justified at the moment; the correct conclusion can be 
made only if the comparison to theory is incorporated into the analysis at early 
stage, and, therefore, it must be left for experiment itself. For this reason in our 
present analysis we aim only at tentatative conclusions and rather try to illustrate 
the sensitivity of the spectra to underlying theoretical parameters; this procedure 
can give an idea of what kind of precision one can expect in determination of these 
parameters from future accurate measurements. In other words, in varying input 
parametes we, to some extent, take more literally the theoretical spectrum obtained 
for the central values of parameters as representing the possible experimental shape, 
and compare other spectra with this hypothetical one. 
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In Fig. 4 we show the full spectrum in its high energy part varying the b quark 
mass: mb(0.7GeV) = 4.67 GeV (solid line), 4.58 GeV (dashed line) and 4.83 GeV 
(long-dashed line); these (rather odd at first sight) values are the same as in Fig. 3 
where they merely correspond to the different choice of the normalization point. 
We stress that now the normalization point is fixed at 0.7 GeV. We choose to take 
the examples above to separately study the effects of independent varying mass 
or normalization point. The value of /x^ is kept fixed at 0.5 GeV^. Experimental 
points are also shown, with arbitrary normalization; the normalization of all three 
theoretical spectra is the same. 

In Fig. 5, on the contrary, we fix m;,(0.7GeV) = 4.67 GeV and try to vary 
between 0.25 GeV^ (dashed line) and 0.75 GeV^ (long-dashed line). 

A look at Figs. 4 and 5 suggests that at present all theoretically allowed range 
of parameters is consistent with the CLEO data |l| , and, apparently, even noticeably 
wider range of parametes cannot be excluded. Still it seems that the most probable 
values of rrib and /i^ are likely to lie not far from our current theoretical expectation. 
Although, clearly, any more definite conclusion can emerge only from a dedicated 
experimental fit. 

In relation to the future data. Fig. 4 shows that the endpoint spectrum is rather 
sensitive to the value of the quark mass, and it seems probable that the precise data 
on the photon spectrum will allow to determine it here with the accuracy of about 
50 MeV. On the other hand, according to Fig. 5 the dependence of the shape on 
the expectation value of the kinetic operator is not too strong and the possibility 
to decrease the error in fx^ below 0.15 GeV^ looks questionable, if one takes into 
account that the theoretical resolution in the photon energy one can pretend on in 
this approach, is not more than ~ 50 MeV. 



7 Conclusions 

Shortly after it was realized that the parton-model delta function in the b —>■ sj 
spectrum is substituted by a universal primordial distrubition in B —>■ Xgj concerns 
were aired that the hard gluon radiation of the Sudakov type might totally wash 
away the endpoint peak, so that no traces of the primordial distribution will be 
visible. The main lesson we learn from the present work is that this is not the 
case. Although the hard gluon emission definitely smears the primordial function, 
when properly treated within OPE this perturbative smearing turns out to be quite 
modest. As is seen from the plots discussed in Sect. 5 the peak in the endpoint 
domain survives perturbative smearing and is quite conspicuous. We hasten to add 
that it is not an artefact of the finite value of rrib used in the numerical analysis. 



The survival of the primordial peak was anticipated in Ref. |jT2 . 

Several reservations must be made as to the accuracy of the various approxima- 
tions we adopted. Being QCD-compatible, the theory we suggest is simplified in 
many aspects. 
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First, in the OPE analysis we dealt with the leading twist operators only. Higher 
twist operators will introduce 1/mi, corrections, generally speaking, in all relations 
derived above, which may or may not be numerically significant. We know for sure, 
however, that the higher twist operators are absolutely crucial in a part of the win- 
dow, namely in the region which will be called resonance, Mb/2 — ~ A^qj^/ irtb. 
By definition, in this domain of the spectrum the hadronic system produced in the 
radiative B decay is just one low-lying strange resonance, like K*. In the discussion 
above we ignored the existence of this region altogether, consistently repeating that 
the physical spectrum stretches up to Mb/2. This was not wrong in our approxi- 
mation considering that the resonance domain is an effect showing up only at the 
(relative) resolution of l/rn^. However, since is not academically large in prac- 
tical terms it is impossible to ignore the existence of the resonance region, at least 
for pure kinematical reasons. The mechanism governing the hadronic transition in 
the resonance domain may be quite different - a form factor type mechanism 



Then our prediction for the physical spectrum can not be valid point-by-point on 
the extreme right of the spectrum, when E approaches the kinematic boundary by 
distance of order K^Q-^/mi,. This resonance domain occupies, roughly, ~ 1/5 of the 
window. 

In general, since only the leading twist operators were considered, in this approx- 
imation one cannot, in principle, predict the spectrum with the energy resolution of 
the order of AQQj-,/mfe. Before confronting theory and experiment a smearing over 
a larger interval of energy must be perfomed. The smearing is superfluous outside 
the resonance domain where the measured spectrum is smooth by itself. On the 
other hand, if there is any hope of extrapolating the prediction for the spectrum 
to the resonance domain, the prediction will refer to the smeared spectrum, not 
point-by-point. 

With very little information about the high twist contributions, one can get an 
idea about the size of the smearing interval resorting to consideration of the impact 
of the p — TT or K* — K mass difference on the kinematics of the decay |[T^. (This 
splitting is the most transparent effect of the next-to-leading twist corrections.) The 
corresponding difference in the two body decay energy 5E is given by 

Ml, - Ml Ml - Ml 

6E ~ ~ ^ ~ 60 MeV . (42) 

2mi, 2m}j 

We see, therefore, that the size of the resonance domain and the smearing interval, 
although formally scaling like AQc^/mfc, is quite signiflcant numerically. 

It may have important practical implications. For example, the point-to-point 
extraction of the distribution function from experimental data in the region very 
close to the kinematic boundary, sometimes proposed for elimination of the model 
dependence of the endpoint determination of Vuh, is not likely to be accurate enough 
for that reason. 

On the other hand, the effect of inclusion of the current quark mass in the 
analysis is minute and can be safely neglected: it scales like m1/2mb (or is linear 
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in rris in higher twists), and is much smaller numerically than the effects discussed 
above. 

Addressing the phenomenological implications of our analysis, we report that in 
spite of the crude nature of the current data for the photon spectrum in B —>■ Xg'j, it 
is still possible to conclude that the current theory including both, perturbative and 
nonperturbative effects, is consistent. Moreover, the existing theoretical estimates 
for nib and fi^ 0, ^ generate the spectrum which seems rather close to the 
data, rather than lying near the edge of contradicting them. Bearing in mind a very 
crude character of the existing data points we did not attempt any fit of the data 
that would allow us to extract mi,, and /i^. With data becoming more accurate such 
a fit will become a necessity. 

We emphasize that our analysis in the future, when more exact data will appear, 
will change only in a technical sense. First, one needs to go beyond the "natural" 
approximation used in the present paper and rely on APS, as it is described in 
Sects. 4.5-4.6. The full 0{as) calculation of the effective 6^5 + 7 vertex is also 
desirable if one wants to consider the whole spectrum and not only the endpoint 
region, for example, address the total decay width on a parallel footing. Second, one 
may need to consider a more general ansatz for the primordial distribution function 
F{x), say, put k = 3, which allows a much wider variety of shapes of the primordial 
distribution. Then even some experimental estimates of the third moment may be 
obtained. 

With new precise experimental data one can hope to get a rather good indepen- 
dent estimate of the running mass at a relatively low normalization point with an 
accuracy as high as ~ 50 MeV. Also, a reasonable estimate will be possible for the 
kinetic matrix element yU^; the accuracy of ~ 0.1 GeV^ is conceivable here. We must 
note that the primordial distribution (depending on the low-energy parameters) is 
still rather significantly smeared by the short- distance perturbative effects. For this 
reason the apparent effect of varying the low-energy parameters is softened and their 
determination would require a dedicated analysis. The intervention of short distance 
corrections is much less pronounced in the semileptonic b civ transitions and, in 
particular, in the SV kinematics. The latter process, thus, offers a more accurate 
and less involoved detemination of rrih and /i^ along the suggestions of Refs. [Bl, P3| . 



The corresponding analysis is under way now ||69|. Consideration of the semileptonic 
h — > ulv transitions is very similar to what is done here for h 57. This work is also 
in progress [^]. It seems important to make independent determinations in different 
situations of heavy-to-light and heavy-to-heavy transitions where the corresponding 
distribution functions are radically different but are related to each other by certain 



constraints 12 



From the theoretical standpoint, in the present paper we develop the OPE- 
based formalism, going explicitly beyond the practical version of OPE, using the 
inclusive 6 ^ 57 transition as an example. For practical purposes we accepted the 
approximation which basically is similar to the double logarithmic approximation 
but incorpoprates the BLM-type improvement accounting for the perturbative run- 
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ning of the strong coupling to all orders; we call it "natural" approximation. Our 
consideration seems instructive in revealing several key elements of the OPE-based 
approach. One clearly sees here the necessity of introduction of the separation scale 
fi to discriminate between short- distance and long-distance effects. The discrimina- 
tion is rather nontrivial in the Minkowski kinematics, even in calculating the purely 
perturbative corrections. On the other hand, it is shown that using the well defined 
short distance parameters like mb{fi) which do depend in the calculable way on fi, 
rather than ill defined "absolute" parameters of the sort of "resummed" pole mass 
originally suggested in HQET, leads to a consistent and //-independent physical 
spectrum. 

We also illustrated that seeds of the basic features of the physical spectrum - 
usually these features are attributed purely to nonperturbative effects - are in fact 
seen in the perturbative calculations as long as the latter are done in the OPE 
consistent manner. In particular, we showed that a "window" which manifests itself 
as the interval of the physical spectrum stretching above the kinematical boundary 
possible for the parton level process, is present already in the proper perturbative 
treatment. The consistent application of the "natural" expressions also allowed 
us to estimate the asymptotic behavior of the perturbative end point spectrum at 
(academically) high mass of the decaying heavy quark, which is /i-dependent as well. 

The more accurate description of the spectrum requires going beyond the dou- 
ble logarithmic approximation, even improved by accounting for the running of the 
strong coupling. In particular, the question of the renormalization point dependence 
of the kinetic operator requires specifying explicitly the normalization scheme; the 
difference between the schemes, however, is beyond the double logarithmic accu- 
racy in the description of the spectrum. Numerically such effects constitute about 
0.1 GeV^ in /i^ [^. We suggest a more advanced and, on the other hand, techni- 



cally rather feasible approximation - the APS. There are good reasons to believe 
that APS must yield sufficient numerical accuracy: it incorporates the exact 0{as) 
result, the summation of the Sudakov logarithms and, additionally, the running of 
as within the BLM hypothesis. The numerical corrections to this approximation, 
therefore, are expected to be rather small for b decays. 

Another new theoretical element of the corresponding analysis is that calcu- 
lation of the perturbative effects is performed in the way required by OPE, viz., 
the perturbative coefficient functions are obtained using the explicit infrared cutoff 
fi. This procedure, more or less straightforward in the Euclidean calculations, is 
rather involved technically in all Minkowski processes where spectral densities at 
nonperturbative level are addressed. No consistent approach existed in the litera- 
ture previously even in the framework of the BLM approximation. This problem 
must be solved in order to construct APS. We showed explicitly how it can be done 
with the example of the particular regularization scheme which is very convenient 
in the BLM-type calculations; no problems with infrared renormalons and related 
ambiguites in resummation of divergent series appear here as a manifestation of 
the consistent application of Wilson's OPE, and nonperturbative matrix elements 
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have definite meaning. This method is apphcable in the direct way to all quantities 
treated within OPE. In the context of the present paper we had to limit ourselves 
only to a brief description of the method and formulating some results; more details 
and the discussion of both physical consequences and technical questions will be 
given elsewhere pU |. 
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Figure Captions 



Fig. 1 The region of the values of 02 and possible with the ansatz for F{x) given 
by Eq. (0) with k = 1 (solid line) and = 3 (dashed line). We plot the value of 
as — a™°|fc=i along the y axis, where a™™|fc=i = ~2a\. 

Fig. 2 The primordial distribution (dashed line), the perturbative spectrum (dotted 
line) and the full spectrum (solid line) obtained for yU = 0.7 GeV, mb(/i) = 4.67 GeV, 
= 0.5 GeV^ and Aq^j-, = 300 MeV. All curves shown have the same area; the 
height of the elastic peak in the perturbative spectrum corresponds to its weight of 
0.454 in the total width. Further details are described in Sect. 6. 

Fig. 3 The full end point spectrum obtained at different values of the normalization 
point fi = 0.7 GeV (sohd line), 1.0 GeV (dashed line) and 0.4 GeV (long-dashed 
line); the running masses mb{fi) are used, mfe(0.7GeV) = 4.67 GeV, mf,(l GeV) = 
4.58 GeV and mfe(0.4 GeV) = 4.83 GeV . The value of fxl = 0.5 GeV^ is kept fixed. 

Fig. 4 The dependence of the end point spectrum on the heavy quark mass; the 
normalization point is /x = 0.7 GeV. Solid line is ruf, = 4.67 GeV, dashed line 
represents rrih = 4.58 GeV and long-dashed line corresponds to ruf, = 4.83 GeV. The 
value of fj,"^ = 0.5 GeV^ is taken for all curves. CLEO data are shown in an 
aribitrary normalization. 

Fig. 5 The dependence of the end point spectrum on the expectation value fi^. of 
the kinetic operator: solid line is /i^ = 0.5 GeV^, dashed line is fx^ = 0.25 GeV^ and 
long-dashed lined is /i^ = 0.75 GeV^. We assume mi, = 4.67 GeV and /i = 0.7 GeV. 
CLEO data [Q] are shown in an aribitrary normalization. 
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